DeTox: a pipeline for the detection of toxins in venomous organisms.

Venomous organisms have independently evolved the ability to produce toxins 101 times during their evolutionary history, resulting in over 200 000 venomous species. Collectively, these species produce millions of toxins, making them a valuable resource for bioprospecting and understanding the evolutionary mechanisms underlying genetic diversification. RNA-seq is the preferred method for characterizing toxin repertoires, but the analysis of the resulting data remains challenging. While early approaches relied on similarity-based mapping to known toxin databases, recent studies have highlighted the importance of structural features for toxin detection. The few existing pipelines lack an integration between these complementary approaches, and tend to be difficult to run for non-experienced users. To address these issues, we developed DeTox, a comprehensive and user-friendly tool for toxin research. It combines fast execution, parallelization and customization of parameters. DeTox was tested on published transcriptomes from gastropod mollusks, cnidarians and snakes, retrieving most putative toxins from the original articles and identifying additional peptides as potential toxins to be confirmed through manual annotation and eventually proteomic analysis. By integrating a structure-based search with similarity-based approaches, DeTox allows the comprehensive characterization of toxin repertoire in poorly-known taxa. The effect of the taxonomic bias in existing databases is minimized in DeTox, as mirrored in the detection of unique and divergent toxins that would have been overlooked by similarity-based methods. DeTox streamlines toxin annotation, providing a valuable tool for efficient identification of venom components that will enhance venom research in neglected taxa.

Taxonomic Distribution and Molecular Evolution of Mytilectins.

R-type lectins are a widespread group of sugar-binding proteins found in nearly all domains of life, characterized by the presence of a carbohydrate-binding domain that adopts a ß-trefoil fold. Mytilectins represent a recently described subgroup of ß-trefoil lectins, which have been functionally characterized in a few mussel species (Mollusca, Bivalvia) and display attractive properties, which may fuel the development of artificial lectins with different biotechnological applications. The detection of different paralogous genes in mussels, together with the description of orthologous sequences in brachiopods, supports the formal description of mytilectins as a gene family. However, to date, an investigation of the taxonomic distribution of these lectins and their molecular diversification and evolution was still lacking. Here, we provide a comprehensive overview of the evolutionary history of mytilectins, revealing an ancient monophyletic evolutionary origin and a very broad but highly discontinuous taxonomic distribution, ranging from heteroscleromorphan sponges to ophiuroid and crinoid echinoderms. Moreover, the overwhelming majority of mytilectins display a chimera-like architecture, which combines the ß-trefoil carbohydrate recognition domain with a C-terminal pore-forming domain, suggesting that the simpler structure of most functionally characterized mytilectins derives from a secondary domain loss.

Deciphering the role of recurrent FAD-dependent enzymes in bacterial phosphonate catabolism.

Phosphonates-compounds containing a direct C-P bond-represent an important source of phosphorus in some environments. The most common natural phosphonate is 2-aminoethylphosphonate (AEP). Many bacteria can break AEP down through specialized "hydrolytic" pathways, which start with the conversion of AEP into phosphonoacetaldehyde (PAA), catalyzed by the transaminase PhnW. However, the substrate scope of these pathways is very narrow, as PhnW cannot process other common AEP-related phosphonates, notably N-methyl AEP (M1AEP). Here, we describe a heterogeneous group of FAD-dependent oxidoreductases that efficiently oxidize M1AEP to directly generate PAA, thus expanding the versatility and usefulness of the hydrolytic AEP degradation pathways. Furthermore, some of these enzymes can also efficiently oxidize plain AEP. By doing so, they surrogate the role of PhnW in organisms that do not possess the transaminase and create novel versions of the AEP degradation pathways in which PAA is generated solely by oxidative deamination.

Relaxation of Natural Selection in the Evolution of the Giant Lungfish Genomes.

Nonadaptive hypotheses on the evolution of eukaryotic genome size predict an expansion when the process of purifying selection becomes weak. Accordingly, species with huge genomes, such as lungfish, are expected to show a genome-wide relaxation signature of selection compared with other organisms. However, few studies have empirically tested this prediction using genomic data in a comparative framework. Here, we show that 1) the newly assembled transcriptome of the Australian lungfish, Neoceratodus forsteri, is characterized by an excess of pervasive transcription, or transcriptional leakage, possibly due to suboptimal transcriptional control, and 2) a significant relaxation signature in coding genes in lungfish species compared with other vertebrates. Based on these observations, we propose that the largest known animal genomes evolved in a nearly neutral scenario where genome expansion is less efficiently constrained.

Gene presence/absence variation in Mytilus galloprovincialis and its implications in gene expression and adaptation.

Presence/absence variation (PAV) is a well-known phenomenon in prokaryotes that was described for the first time in bivalves in 2020 in Mytilus galloprovincialis. The objective of the present study was to further our understanding of the PAV phenomenon in mussel biology. The distribution of PAV was studied in a mussel chromosome-level genome assembly, revealing a widespread distribution but with hotspots of dispensability. Special attention was given to the effect of PAV in gene expression, since dispensable genes were found to be inherently subject to distortions due to their sparse distribution among individuals. Furthermore, the high expression and strong tissue specificity of some dispensable genes, such as myticins, strongly supported their biological relevance. The significant differences in the repertoire of dispensable genes associated with two geographically distinct populations suggest that PAV is involved in local adaptation. Overall, the PAV phenomenon would provide a key selective advantage at the population level.

Biochemical and Functional Characterization of the Three Zebrafish Transglutaminases 2.

Transglutaminase 2 (TG2) is a multifunctional protein widely distributed in various tissues and involved in many physiological and pathological processes. However, its actual role in biological processes is often controversial as TG2 shows different effects in these processes depending on its localization, cell type, or experimental conditions. We characterized the enzymatic and functional properties of TG2 proteins expressed in Danio rerio (zebrafish) to provide the basis for using this established animal model as a reliable tool to characterize TG2 functions in vivo. We confirmed the existence of three genes orthologous to human TG2 (zTGs2) in the zebrafish genome and their expression and function during embryonic development. We produced and purified the zTGs2s as recombinant proteins and showed that, like the human enzyme, zTGs2 catalyzes a Ca2+ dependent transamidation reaction that can be inhibited with TG2-specific inhibitors. In a cell model of human fibroblasts, we also demonstrated that zTGs2 can mediate RGD-independent cell adhesion in the extracellular environment. Finally, we transfected and selected zTGs2-overexpressing HEK293 cells and demonstrated that intracellular zTGs2 plays a very comparable protective/damaging role in the apoptotic process, as hTG2. Overall, our results suggest that zTGs2 proteins behave very similarly to the human ortholog and pave the way for future in vivo studies of TG2 functions in zebrafish.

Bivalves Present the Largest and Most Diversified Repertoire of Toll-Like Receptors in the Animal Kingdom, Suggesting Broad-Spectrum Pathogen Recognition in Marine Waters.

Toll-like receptors (TLRs) are the most widespread class of membrane-bound innate immune receptors, responsible of specific pathogen recognition and production of immune effectors through the activation of intracellular signaling cascades. The repertoire of TLRs was analyzed in 85 metazoans, enriched on molluscan species, an underrepresented phylum in previous studies. Following an ancient evolutionary origin, suggested by the presence of TLR genes in Anthozoa (Cnidaria), these receptors underwent multiple independent gene family expansions, the most significant of which occurred in bivalve molluscs. Marine mussels (Mytilus spp.) had the largest TLR repertoire in the animal kingdom, with evidence of several lineage-specific expanded TLR subfamilies with different degrees of orthology conservation within bivalves. Phylogenetic analyses revealed that bivalve TLR repertoires were more diversified than their counterparts in deuterostomes or ecdysozoans. The complex evolutionary history of TLRs, characterized by lineage-specific expansions and losses, along with episodic positive selection acting on the extracellular recognition domains, suggests that functional diversification might be a leading evolutionary force. We analyzed a comprehensive transcriptomic data set from Mytilus galloprovincialis and built transcriptomic correlation clusters with the TLRs expressed in gills and in hemocytes. The implication of specific TLRs in different immune pathways was evidenced, as well as their specific modulation in response to different biotic and abiotic stimuli. We propose that, in a similar fashion to the remarkable functional specialization of vertebrate TLRs, the expansion of the TLR gene family in bivalves attends to a functional specification motivated by the biological particularities of these organisms and their living environment.

Gene expression profiling in white blood cells reveals new insights into the molecular mechanisms of thalidomide in children with inflammatory bowel disease.

Thalidomide has emerged as an effective immunomodulator in the treatment of pediatric patients with inflammatory bowel disease (IBD) refractory to standard therapies. Cereblon (CRBN), a component of E3 protein ligase complex that mediates ubiquitination and proteasomal degradation of target proteins, has been identified as the primary target of thalidomide. CRBN plays a crucial role in thalidomide teratogenicity, however it is unclear whether it is also involved in the therapeutic effects in IBD patients. This study aimed at identifying the molecular mechanisms underpinning thalidomide action in pediatric IBD. In this study, ten IBD pediatric patients responsive to thalidomide were prospectively enrolled. RNA-sequencing (RNA-seq) analysis and functional enrichment analysis were carried out on peripheral blood mononuclear cells (PBMC) obtained before and after twelve weeks of treatment with thalidomide. RNA-seq analysis revealed 378 differentially expressed genes before and after treatment with thalidomide. The most deregulated pathways were cytosolic calcium ion concentration, cAMP-mediated signaling, eicosanoid signaling and inhibition of matrix metalloproteinases. Neuronal signaling mechanisms such as CREB signaling in neurons and axonal guidance signaling also emerged. Connectivity Map analysis revealed that thalidomide gene expression changes were similar to those exposed to MLN4924, an inhibitor of NEDD8 activating enzyme, suggesting that thalidomide exerts its immunomodulatory effects by acting on the ubiquitin-proteasome pathway. In vitro experiments on cell lines confirmed the effect of thalidomide on candidate altered pathways observed in patients. These results represent a unique resource for enhanced understanding of thalidomide mechanism in pediatric patients with IBD, providing novel potential targets associated with drug response.

Unveiling the Impact of Gene Presence/Absence Variation in Driving Inter-Individual Sequence Diversity within the CRP-I Gene Family in Mytilus spp.

Mussels (Mytilus spp.) tolerate infections much better than other species living in the same marine coastal environment thanks to a highly efficient innate immune system, which exploits a remarkable diversification of effector molecules involved in mucosal and humoral responses. Among these, antimicrobial peptides (AMPs) are subjected to massive gene presence/absence variation (PAV), endowing each individual with a potentially unique repertoire of defense molecules. The unavailability of a chromosome-scale assembly has so far prevented a comprehensive evaluation of the genomic arrangement of AMP-encoding loci, preventing an accurate ascertainment of the orthology/paralogy relationships among sequence variants. Here, we characterized the CRP-I gene cluster in the blue mussel Mytilus edulis, which includes about 50 paralogous genes and pseudogenes, mostly packed in a small genomic region within chromosome 5. We further reported the occurrence of widespread PAV within this family in the Mytilus species complex and provided evidence that CRP-I peptides likely adopt a knottin fold. We functionally characterized the synthetic peptide sCRP-I H1, assessing the presence of biological activities consistent with other knottins, revealing that mussel CRP-I peptides are unlikely to act as antimicrobial agents or protease inhibitors, even though they may be used as defense molecules against infections from eukaryotic parasites.

Chromosome-Level Genome Assembly of the Blue Mussel Mytilus chilensis Reveals Molecular Signatures Facing the Marine Environment.

The blue mussel Mytilus chilensis is an endemic and key socioeconomic species inhabiting the southern coast of Chile. This bivalve species supports a booming aquaculture industry, which entirely relies on artificially collected seeds from natural beds that are translocated to diverse physical-chemical ocean farming conditions. Furthermore, mussel production is threatened by a broad range of microorganisms, pollution, and environmental stressors that eventually impact its survival and growth. Herein, understanding the genomic basis of the local adaption is pivotal to developing sustainable shellfish aquaculture. We present a high-quality reference genome of M. chilensis, which is the first chromosome-level genome for a Mytilidae member in South America. The assembled genome size was 1.93 Gb, with a contig N50 of 134 Mb. Through Hi-C proximity ligation, 11,868 contigs were clustered, ordered, and assembled into 14 chromosomes in congruence with the karyological evidence. The M. chilensis genome comprises 34,530 genes and 4795 non-coding RNAs. A total of 57% of the genome contains repetitive sequences with predominancy of LTR-retrotransposons and unknown elements. Comparative genome analysis of M. chilensis and M. coruscus was conducted, revealing genic rearrangements distributed into the whole genome. Notably, transposable Steamer-like elements associated with horizontal transmissible cancer were explored in reference genomes, suggesting putative relationships at the chromosome level in Bivalvia. Genome expression analysis was also conducted, showing putative genomic differences between two ecologically different mussel populations. The evidence suggests that local genome adaptation and physiological plasticity can be analyzed to develop sustainable mussel production. The genome of M. chilensis provides pivotal molecular knowledge for the Mytilus complex.

Helminthic host defense peptides: using the parasite to defend the host.

Parasitic helminths are destined to share niches with a variety of microbiota that inevitably influence their interaction with the host. To modulate the microbiome for their benefit and defend against pathogenic isolates, helminths have developed host defense peptides (HDPs) and proteins as integral elements of their immunity. These often exert a relatively nonspecific membranolytic activity toward bacteria, sometimes with limited or no toxicity toward host cells. With a few exceptions, such as nematode cecropin-like peptides and antibacterial factors (ABFs), helminthic HDPs are largely underexplored. This review scrutinizes current knowledge on the repertoire of such peptides in helminths and promotes their research as potential leads for an anti-infective solution to the burgeoning problem of antibiotic resistance.

A high-quality reference genome for the critically endangered Aeolian wall lizard, Podarcis raffonei.

The Aeolian wall lizard, Podarcis raffonei, is an endangered species endemic to the Aeolian archipelago, Italy, where it is present only in 3 tiny islets and a narrow promontory of a larger island. Because of the extremely limited area of occupancy, severe population fragmentation and observed decline, it has been classified as Critically Endangered by the International Union for the Conservation of Nature (IUCN). Using Pacific Biosciences (PacBio) High Fidelity (HiFi) long-read sequencing, Bionano optical mapping and Arima chromatin conformation capture sequencing (Hi-C), we produced a high-quality, chromosome-scale reference genome for the Aeolian wall lizard, including Z and W sexual chromosomes. The final assembly spans 1.51 Gb across 28 scaffolds with a contig N50 of 61.4 Mb, a scaffold N50 of 93.6 Mb, and a BUSCO completeness score of 97.3%. This genome constitutes a valuable resource for the species to guide potential conservation efforts and more generally for the squamate reptiles that are underrepresented in terms of available high-quality genomic resources.

Coming New Age of Marine Glycomics: The Fundamental, Medical, and Ecological Aspects.

This Special Issue "Marine Glycomics" (https://www.mdpi.com/journal/marinedrugs/special_issues/Marine_Glycomics, accessed on 12 September 2022) provided new approachesand information on bioactive compounds, such as glycans and lectins from marine animals,seaweeds, and microorganisms for the application of clinical therapy and elucidation of thephysiological functions of marine organisms [...].

Expansion and Neofunctionalization of Actinoporin-like Genes in Mediterranean Mussel (Mytilus galloprovincialis).

Pore-forming toxins are an important component of the venom of many animals. Actinoporins are potent cytolysins that were first detected in the venom of sea anemones; however, they are occasionally found in animals other than cnidarians and are expanded in a few predatory gastropods. Here, we report the presence of 27 unique actinoporin-like genes with monophyletic origin in Mytilus galloprovincialis, which we have termed mytiporins. These mytiporins exhibited a remarkable level of molecular diversity and gene presence-absence variation, which warranted further studies aimed at elucidating their functional role. We structurally and functionally characterized mytiporin-1 and found significant differences from the archetypal actinoporin fragaceatoxin C. Mytiporin-1 showed weaker permeabilization activity, no specificity towards sphingomyelin, and weak activity in model lipid systems with negatively charged lipids. In contrast to fragaceatoxin C, which forms octameric pores, functional mytiporin-1 pores on negatively charged lipid membranes were hexameric. Similar hexameric pores were observed for coluporin-26 from Cumia reticulata and a conoporin from Conus andremenezi. This indicates that also other molluscan actinoporin-like proteins differ from fragaceatoxin C. Although the functional role of mytiporins in the context of molluscan physiology remains to be elucidated, the lineage-specific gene family expansion event that characterizes mytiporins indicates that strong selective forces acted on their molecular diversification. Given the tissue distribution of mytiporins, this process may have broadened the taxonomic breadth of their biological targets, which would have important implications for digestive processes or mucosal immunity.

Transcriptome Analysis Reveals Early Hemocyte Responses upon In Vivo Stimulation with LPS in the Stick Insect Bacillus rossius (Rossi, 1788).

Despite a growing number of non-model insect species is being investigated in recent years, a greater understanding of their physiology is prevented by the lack of genomic resources. This is the case of the common European stick insect Bacillus rossius (Rossi, 1788): in this species, some knowledge is available on hemocyte-related defenses, but little is known about the physiological changes occurring in response to natural or experimental challenges. Here, the transcriptional signatures of adult B. rossius hemocytes were investigated after a short-term (2 h) LPS stimulation in vivo: a total of 2191 differentially expressed genes, mostly involved in proteolysis and carbohydrate and lipid metabolic processes, were identified in the de novo assembled transcriptome and in-depth discussed. Overall, the significant modulation of immune signals-such as C-type lectins, ML domain-containing proteins, serpins, as well as Toll signaling-related molecules-provide novel information on the early progression of LPS-induced responses in B. rossius.

Independent acquisition of short insertions at the RIR1 site in the spike N-terminal domain of the SARS-CoV-2 BA.2 lineage.

Although the major SARS-CoV-2 omicron lineages share over 30 non-synonymous substitutions in the spike glycoprotein, they show several unique mutations that were acquired after their ancestral split. One of the most intriguing mutations associated with BA.1 is the presence of the inserted tripeptide Glu-Pro-Glu within the N-terminal domain, at a site that had previously independently acquired short insertions in several other SARS-CoV-2 lineages. Although the functional implications of the small nucleotide sequences found at this insertion hotspot, named RIR1, are still unclear, we have previously hypothesized that they may play a compensatory role in counterbalancing minor fitness deficits associated with other co-occurring spike non-synonymous mutations. Here, we show that similar insertion events have independently occurred at RIR1 at least 20 times in early 2022 within the BA.2 lineage, being occasionally associated with significant community transmission. One of these omicron sublineages, characterized by a Ser-Gly-Arg insertion in position 212, has been responsible for over 4000 documented COVID-19 cases worldwide between January and July 2022, for the most part concentrated in Denmark, where it reached a national prevalence close to 4% (10% in the Nordjylland region) in mid-May. Although the concurrent spread of the BA.2.12.1, BA.4 and BA.5 lineages led to the rapid decline of this BA.2 sublineage, the independent acquisition of several other RIR1 insertions on a BA.2 genomic background suggests that these events may provide a slight fitness advantage. Therefore, they should be carefully monitored in the upcoming months in other emerging omicron-related lineages, including BA.5.

Detecting Structural Variants and Associated Gene Presence-Absence Variation Phenomena in the Genomes of Marine Organisms.

As complete genomes become easier to attain, even from previously difficult-to-sequence species, and as genomic resequencing becomes more routine, it is becoming obvious that genomic structural variation is more widespread than originally thought and plays an important role in maintaining genetic variation in populations. Structural variants (SVs) and associated gene presence-absence variation (PAV) can be important players in local adaptation, allowing the maintenance of genetic variation and taking part in other evolutionarily relevant phenomena. While recent studies have highlighted the importance of structural variation in Mollusca, the prevalence of this phenomenon in the broader context of marine organisms remains to be fully investigated.Here, we describe a straightforward and broadly applicable method for the identification of SVs in fully assembled diploid genomes, leveraging the same reads used for assembly. We also explain a gene PAV analysis protocol, which could be broadly applied to any species with a fully sequenced reference genome available. Although the strength of these approaches have been tested and proven in marine invertebrates, which tend to have high levels of heterozygosity, possibly due to their lifestyle traits, they are also applicable to other species across the tree of life, providing a ready means to begin investigations into this potentially widespread phenomena.

Comparative analysis of novel and common reference genes in adult tissues of the mussel Mytilus galloprovincialis.

BACKGROUND: Real-time quantitative PCR is a widely used method for gene expression analyses in various organisms. Its accuracy mainly relies on the correct selection of reference genes. Any experimental plan involving real-time PCR needs to evaluate the characteristics of the samples to be examined and the relative stability of reference genes. Most studies in mollusks rely on reference genes commonly used in vertebrates. RESULTS: In this study, we focused on the transcriptome of the bivalve mollusk Mytilus galloprovincialis in physiological state to identify suitable reference genes in several adult tissues. Candidate genes with highly stable expression across 51 RNA-seq datasets from multiple tissues were selected through genome-wide bioinformatics analysis. This approach led to the identification of three genes (Rpl14, Rpl32 and Rpl34), whose suitability was evaluated together with 7 other reference genes commonly reported in literature (Act, Cyp-A, Ef1α, Gapdh, 18S, 28S and Rps4). The stability analyses performed with geNorm, NormFinder and Bestkeeper identified specific either single or pairs of genes suitable as references for gene expression analyses in specific tissues and revealed the Act/Cyp-A pair as the most appropriate to analyze gene expression across different tissues. CONCLUSION: Mytilus galloprovincialis is a model system increasingly used in ecotoxicology and molecular studies. Our transcriptome-wide approach represents the first comprehensive investigation aimed at the identification of suitable reference genes for expression studies in this species.

The Mantle Transcriptome of Chamelea gallina (Mollusca: Bivalvia) and Shell Biomineralization.

The striped venus Chamelea gallina is a bivalve mollusc that represents one of the most important fishery resources of the Adriatic Sea. In this work, we investigated for the first time the ability of this species to modulate the expression of genes encoding proteins involved in biomineralization process in response to biotic and abiotic factors. We provided the first comprehensive transcriptome from the mantle tissue of clams collected in two sampling sites located along the Italian Adriatic coast and characterized by different environmental features. Moreover, the assessment of environmental parameters, scanning electron microscopy (SEM), and X-ray diffraction (XRD) measurements on valves were conducted to better contextualize RNA sequencing (RNA-Seq) data. Functional annotation of differentially expressed genes (DEGs) and SEM observations highlighted a different shell mineralization behaviour in C. gallina clams collected from two selected sites characterized by diverse environmental parameters.

Transcriptional Contribution of Transposable Elements in Relation to Salinity Conditions in Teleosts and Silencing Mechanisms Involved.

Fish are an interesting taxon comprising species adapted to a wide range of environments. In this work, we analyzed the transcriptional contribution of transposable elements (TEs) in the gill transcriptomes of three fish species exposed to different salinity conditions. We considered the giant marbled eel Anguilla marmorata and the chum salmon Oncorhynchus keta, both diadromous, and the marine medaka Oryzias melastigma, an euryhaline organism sensu stricto. Our analyses revealed an interesting activity of TEs in the case of juvenile eels, commonly adapted to salty water, when exposed to brackish and freshwater conditions. Moreover, the expression assessment of genes involved in TE silencing mechanisms (six in heterochromatin formation, fourteen known to be part of the nucleosome remodeling deacetylase (NuRD) complex, and four of the Argonaute subfamily) unveiled that they are active. Finally, our results evidenced for the first time a krüppel-associated box (KRAB)-like domain specific to actinopterygians that, together with TRIM33, might allow the functioning of NuRD complex also in fish species. The possible interaction between these two proteins was supported by structural prediction analyses.